Role of the Dzyaloshinskii-Moriya interaction in multiferroic perovskites 



I. A. Sergienko and E. Dagotto 
Condensed Matter Sciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA and 
Department of Physics, The University of Tennessee, Knoxville, TN 37996, USA 

With the perovskite multiferroic _RMn03 (R = Gd, Tb, Dy) as guidance, we argue that the 
Dzyaloshinskii-Moriya interaction (DMI) provides the microscopic mechanism for the coexistence 
and strong coupling between ferroelectricity and incommensurate magnetism. We use Monte-Carlo 
simulations and zero temperature exact calculations to study a model incorporating the double- 
exchange, superexchange, Jahn- Teller and DMI terms. The phase diagram contains a multiferroic 
phase between A and E antiferromagnetic phases, in excellent agreement with experiments. 
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I. INTRODUCTION 

Several recent discoveries of unusually strong coupling 
between the ferroelectric (FE) and magnetic order pa- 
rameters have revived the interest in the magnetoelectric 
effect. 1 Due to the possibility of easily controlling the 
electric properties using magnetic fields, a class of com- 
pounds, in which the magnetic order is incommensurate 
(1C) with the lattice period, is particularly interesting for 
future applications^^ Surprisingly, this class of multifer- 
roic materials includes compounds with very diverse crys- 
tallographic structures: the perovskite i?MnC>3 (R = Gd, 
Tb, Dy), 2 i 4 i 5 i 6 i 7 i 8 orthorhombic i?Mn 2 5 (R = Tb, Ho, 
Dy) ) 9 ! 10 !' 1 ! 12 ! 13 hexagonal Ba .5Sri.5Zn 2 Fei 2 O 22f ±i and 
Kagome-staircase M3V2O8. 15 Unfortunately, the values 
of the electric polarization achieved so far are two or- 
ders of magnitude smaller than those in traditional fer- 
roelectrics. However, it is essential to theoretically under- 
stand the new mechanism of magnetoelectric coupling. 

Generally, certain types of magnetic order can lower 
the symmetry of the system to one of the polar 
groups, which allows for ferroelectricity. This is 
a central argument of the existing phenomenologi- 
cal modelsii 5 ^ 6 - According to the recent experimen- 
tal results in Bao. 5 Sri. 5 Zn 2 Fei 2 22f i4 Ni 3 V 2 8f i 5 . and 
TbMn03, 7 helical magnetic structures are the most likely 
candidates to host ferroelectricity In addition, X-ray 
diffraction studies in a number of the above materials 
have revealed that the modulated magnetic structure is 
accompanied by structural modulation ! 2 ! 8 ! 11 It is, there- 
fore, a natural assumption that lattice displacements ac- 
tively participate in the formation of the ferroelectric 
state as well, even though, the FE displacements have not 
been measured directly, owing to their smallness (~ 10 -3 
A, as can be deduced from the value of the FE polariza- 
tion P). This calls for theoretical microscopic models 
providing a mechanism by which the FE lattice displace- 
ments are induced and coupled to the IC magnetic struc- 
ture. 

In this paper, we concentrate on the perovskite man- 
ganites since their experimental phase diagram has 
been studied in much detail&^£ We show that the 
Dzyaloshinksii-Moriya interaction (DMI), linearly depen- 




FIG. 1: (Color online) (a, b) Sketch of the ground state struc- 
ture of TbMnOs. (a) Projection of the Mn spins onto the xy- 
plane. Dashed lines are the boundaries of the unit cell in the 
orthorhombic setting. The diagram on the right illustrates the 
emergence of the in-plane component of P in applied magnetic 
field, (b) The zz-projection of the spin structure and posi- 
tions occupied by the O ions (filled circles), (c) The Mn2 
"molecule" in the cubic perovskite structure. The vector r n 
denotes the displacement of the O ion. 



dent on the displacements of the oxygen ions surround- 
ing transition metal ions, is an essential ingredient in the 
theory of the magnetoelectric effect in IC magnets. In in- 
dependent workji£ it is suggested that the DMI induces 
the polarization of the electronic orbitals, without the 
involvement of the lattice degrees of freedom. In our al- 
ternative scenario, the effect of the DMI is twofold: it 
induces the FE lattice displacements and helps to stabi- 
lize helical magnetic structures at low temperature. 



II. SEMI-EMPIRICAL DISCUSSION 

The cubic perovskite structure of i?MnC>3 is or- 
thorhombically distorted by the GdFe03-type coopera- 
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tive rotation of the MnOg octahedra. Here we use x, 
y, and z to denote the directions in the pseudo-cubic 
setting, and a, b, and c are used for the directions in 
the orthorhombic structure. The c-axis is parallel to z, 
and the relation between the in-plane axes is shown in 
Fig. HJa). In the following, we denote by ICM-FE the 
phase in which the IC magnetic and FE order coexist. 
In i?Mn0 3 (R = Tb, Dy£2& and, possibly^ Gd), the 
ICM-FE phase is the ground state. The magnetic struc- 
ture has the modulation vector k along the 6-axis within 
the xy-plane, and the planes are antiferromagnetically 
stacked along the z-axis. Therefore, for a given Mn ion, 
the two nearest neighbors in the positive directions of the 
x and y axes, have exactly the same spin, as shown in 
Fig-EIa). We use the index n to enumerate the Mn ions 
along a chain in the x-direction [Fig. ^b)]. According 
to a neutron diffraction study of TbMnOs^ the classical 
low-temperature spin structure of the Mn sublattice can 
be described as^ 



S l n = Sq cos(nO + a>i), 



(1) 



where i = x,y, z, Sq = Sq 
The experiment was insensitive to the difference between 
a x and a z , but it was found that a x = a y . This structure 
is stable for temperatures T < 28 K. The lattice structure 
is also IC with the modulation vector k/ = 2 k>£ The FE 
polarization P is oriented along the z-axis. For 28 K < 
T < 41 K, another IC, collinear magnetic structure is 
found with no FE moment 

Considering the possible modifications of spin inter- 
actions by lattice displacements, we first examine the 
isotropic superexchange. To address ferroelectricity, we 
assume that the positions of the Mn ions are fixed, while 
the O ions can displace from their positions in the para- 
electric phase. We consider a chain of Mn ions in the 
x-direction, and denote v n — {x n ,y n: z n ) the displace- 
ment of the O ion located between the Mn spins S n and 
S, i+ i from its position in the idealized cubic perovskitc 
structure. In that structure, the Mn2 "molecule" shown 
in Fig. ^c) has inversion symmetry. Since the symmet- 
ric superexchange (S„ ■ S„ + i) is invariant under inversion 
(i. e. mere interchange of S„ and S n +i), this coupling can 
only depend on even powers of r„ . The axial symmetry 
of the Mn20 "molecule" implies 
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Hrs^r — — 



Jo + ~J'\\xl + \j'Avl + z l) 



(S n ■ S n +i), 
(2) 

where Jo, J|j, and J' ± are constants, and the higher-order 
terms in r„ are neglected. 

In the actual orthorhombically distorted structure, the 
oxygen ions are displaced, so that 



r„ = (-l)"r + 5r n , 



(3) 



where r is constant (ro is a fraction of lA)^ and <5r„ 
is the additional displacement, associated with the IC 
structure (Sr n are of the order of 10 -3 A, as mentioned 



above). It can now be shown that the portion of the 
Hamiltonian depending on 5r n is 



5H CX = Y,(- 1 ) n+1 [ J '\\ x oSx n + J' ± (y a Sy n 

n 

+z Sz n )](S n ■ S„ + i) + H ch 
where the elastic energy 



(4) 



(5) 



is assumed isotropic for clarity, the stiffness k > 0, and 
the second order terms in Sr n are neglected in the mag- 
netoelastic term. Minimizing SH CX with respect to the 
displacements and using we obtain 

Sz n - (-1)™^ Y So 2 {cos# + cos[(2n + 1)9 + 2a,]}, 

i 

(6) 

and similar expressions for Sx n and 8y n . Hence, this 
model reproduces the observed structural modulation 
with the wave vector 2 k. 2 However, it is insufficient to 
explain the net FE polarization, since ^2 n r„ = exactly. 

On the contrary, the DMI, i. e. anisotropic exchange 
interaction^ 1 ^ S„ x S n+ i changes its sign under in- 
version. Using group theory, we find that the follow- 
ing expression is invariant under all symmetry opera- 
tions of the Mn02 "molecule" in the perovskite structure: 
D a (r„) • [S n x S»+i], where 

D x (r„) = 7(0, -z n ,y n ), D y (r„) = j(z n> 0, -x n ) (7) 

for the Mn-O-Mn bonds along the x and y axes, respec- 
tively. The form of Q can also be obtained^ 3 , by pertur- 
bative calculations within the Anderson-Moriya theory 
of superexchange^ For the Mn chain in the x-direction, 
the portion of the Hamiltonian depending on Sr n , 



<S#DM = ]TD x («Sr„).[S„xS n+1 ] 



is minimized by 



5z n = — SqSq sin#sin(o! x 



(8) 



(9) 



and Sx n = Sy n = 0. The same result is obtained for the 
Mn-chain in the y-dircction. Hence, the displacements 
in this model do not depend on n, leading to a net FE 
polarization along the z axis. Moreover, Eq. (JSJ explains 
why the IC phase with collinear spin structure is para- 
electric. Clearly, the equilibrium value of 5r n vanishes in 
this case. 

In applied magnetic fields, the magnetic structure 
changes through various phase transitions.-^ Thus far, the 
magnetic structure was only determined experimentally 
for the magnetic field of 4 T oriented along the a axis»£ 
It is similar to the zero-field structure, and it does not 
induce the in-plane component of PA However, as can 
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be seen from Eqs. J7J and iJHJ, the in-plane component of 
P can be induced, in general. It is plausible to assume 
that the Mn sites that have the same spin at zero field, 
also have the same spin in the applied uniform fields 
[Fig.^a)]. Hence, the x-component of P is always equal 
in absolute value but opposite in sign to its y-component, 
and the in-plane polarization is always directed along the 
a axisj^ in agreement with the experiments^^ 

Equation 10 allows us to estimate the value of 7. For 
typical phonon frequencies n ss leV/A 2 , and we obtain 
7 ps lmeV/A. Correspondingly, the magnitude of D(r„) 
is of the order of jz n w 0.1 mcV « 1 K, which agrees well 
with the previous experimental estimations for perovskite 
manganites i 2 ^i 2 9i 2 7 



III. STABILIZATION OF THE MAGNETIC 
HELIX 



To further understand the role of the DMI in the sta- 
bilization of the ICM-FE phase, we consider the Hamil- 
tonian 

H = - taA*A+*P° - JnJ2 Si ■ Si 

iaa/3er i 

+ JAF S i ■ S i+a + ° a ( r ) ' t Si X Si +J 
ia ia 

+y £(<& + Q 2 yi ) + H, T + f E E <&* 
i i m 

(10) 

based on the orbitally-degenerate double exchange 
model, 2 - for a two-dimensional square lattice, repre- 
senting a Mn02 layer, with periodic boundary condi- 
tions and one e g electron per Mn 3+ ion. Here d\ is 
the creation operator for the electron on site i, orbital 
-y 2 (a), 3z 2 — r 2 (b) and carrying spin a. The hop- 
ping integrals are given by i* a = — \/3t* b = — v3ib a = 

t. Here- 



3f 



hi, 



t and ty = V3t* = V3tl 



ab 

ot bb 



after, t is taken as the energy unit. Jh is the Hund's 
coupling constant between the e g electrons with spin 

s i = Z)ck7<7' dJ laa cr aa fd iaa > (<t = Pauli matrices) and 
three ti g electrons treated as a classical three-dimensional 
spin Si. In the following, Jh is assumed to be infinite. 
Jaf > is the isotropic superexchange constant between 
the t2 g spins. For simplicity, we consider DMI only be- 
tween the t2 g spins and assume that O ions can only move 
within the plane. Q m i denote the classical phonon coor- 
dinates associated with the displacements of the four O 
atoms surrounding the Mn site i. The doubly degenerate 
FE mode (Q x i,Q y i) is shown in Fig. [21 The Jahn- Teller 
(JT) interaction term is defined as usual, 2 ^' 29 



H]T = X(qiiPi + q2\T x \ + l?3iTzi) 



2^ 



(2^ + 4 + 4), 
(11) 

where Pi = £ Q(T d\ aa d iarn r ; = £ Q/3(T 4*<r CT a/3 d i0a is the 
orbital pseudospin. The three-dimensional JT phonon 
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FIG. 2: (Color online) Displacements of oxygen ions (large 
circles) corresponding to the FE phonon modes in two dimen- 
sions. 
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FIG. 3: (Color online) Low temperature results for 
model llUt with m = 1, K2 = 10. (a) Ground state phase 
diagram. Broken lines are the MC results for the 8x8 cluster. 
Solid lines are the results of the calculations in the thermo- 
dynamic limit, (b) Typical MC electronic density of states 
(DOS) as a function of the electronic energy uj in the ICM- 
FE phase, (c) Typical MC snapshot in the ICM-FE phase. 
Shown are the ICM modulation vector k and polarization P. 



modes </i,2,3i have been appropriately redefined for our 
model in terms of the two dimensional phonons Q m \- 

The remnant phonon modes (not FE or JT) are as- 
sumed to have the same spring constant «2, an assump- 
tion that can be easily removed if necessary. If k\ = «2, 
the ICM-FE phase is degenerate with other spin-canted 
structures. We assume k\ < K2 in order to lift the de- 
generacy in favor of the ICM-FE phase. 

We have performed extensive Monte-Carlo (MC) sim- 
ulations of H with A = on an 8 x 8 cluster, which is 
sufficiently large for our purposes. The low-temperature 
(T = 0.01) results are presented in Fig. |3| For small 
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7, the Jaf _ 7 phase diagram incorporates the FM phase 
(which is the two-dimensional precursor of the three- 
dimensional A-type antiferromagnetic phase) and the E- 
type phase, in agreement with previous results for smaller 
clusters m& We find that the FM phase is represented by 
two states, separated by the dotted line in Fig.OJa). The 
regular FM phase with all spins aligned in the same di- 
rection is stable on the left hand side of the dotted line, 
while the right hand side is the domain of stability of the 
so-called "twisted" phase with the period equal to the 
cluster length. This is a well-known effect of the finite 
cluster size and periodic boundary conditionsi22*2ii2£ 

Our most important numerical result is the stabiliza- 
tion of a new modulated (ICM-FE) phase at finite 7. The 
presence of the three phases in the ground state phase di- 
agram is in excellent agreement with the experiments 
The period of the ICM-FE phase is fixed in the MC sim- 
ulations due to the small cluster, but it is expected to 
vary smoothly with the parameters of H when the lattice 
size is increased, as discussed below. A typical snapshot 
showing the spin structure and O displacements in this 
phase is presented in Fig. Etc) showing the spin modula- 
tion with k = (7r/2,7r/2). In the two-dimensional model 
studied here, P as well as the spins lie within the plane. 
This is a result of confining the O ions to move within 
the plane. The ground state with P along the z-direction 
can be obtained if model (|10[) is considered in three di- 
mensions and the parameters of the model are changed 
in accordance with the orthorhombic distortions of the 
cubic lattice to remove the equivalence of the crystallo- 
graphic directions. 

Figure E^b) shows the calculated electronic density of 
states in the ICM-FE phase, with a gap at the chemical 
potential /1. This demonstrates the insulating character 
of the new phase, which insures that the polarization is 
not screened out by the free charge carriers. 

To elucidate the influence of the finite size effects on 
the ground state phase diagram, we calculated the ener- 
gies for the idealized spin and phonon configurations, for 
which the e g part of H can be diagonalized exactly in the 
thermodynamic limit. This calculation becomes exact at 
zero temperature for the model considered here, with adi- 
abatic phonons and classical spins, under the reasonable 
assumption that no other phases appear as ground states. 
In particular, we obtain for the ICM-FE phase, 

9 -y 2 

-^icm-fe = N(- 1.60 cos - + 2J AF cos6 — sin 2 6>), 

2 4k 1 

(12) 

where N is the number of Mn sites in the lattice and 9 is 
the angle between two neighboring spins. As can be seen 
from Fig. [21a), while the boundary between the FM and 
E phases remains practically unchanged, the increase of 
the lattice size strongly affects the ICM-FE phase. This 
is expected, since 9 can now assume the continuous val- 
ues found by the minimization of i|12|) , instead of being 
restricted by the boundary conditions. Generally, the ab- 
solute value of k changes from zero at the FM-ICM phase 
transition to some finite value at the discontinuous ICM- 
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FIG. 4: Results of exact diagnolization of Hamiltonian HUH 
at T=0 in the thermodynamic limit for A = 2.0, jci = 1, and 
Jt2 = 10. (a) Ground state phase diagram, (b) Electronic 
density of states in the ICM-FE phase. 

E phase transition. Also, in the phase diagram obtained 
by this procedure, only the regular FM paraelectric phase 
is found for small 7 and Jaf < 0.18. 

The minimal value of 7 needed to stabilize the ICM-FE 
phase between A and E phases is calculated in the ther- 
modynamic limit to be about 0.3 [see Fig. Eta)], which 
corresponds to approximately 200 meV/A in the physical 
units. This is two orders of magnitude larger than the ex- 
pected empirical value obtained in Sec. II. However, the 
numbers can change as the model is improved. In partic- 
ular, our zero temperature results demonstrate that the 
minimal 7 needed can be made infinitesimal if the JT 
interaction is taken into account. In Fig. 0] we show the 
results for a realistic value of A = 2.0. The JT interac- 
tion also significantly improves insulating properties of 
the model. The details of this study for finite A, includ- 
ing the results on orbital order and MC simulations, will 
be reported elsewhere. 

The results of the finite temperature MC study are 
reported in Fig. El Figure EJa) depicts a typical temper- 
ature dependence of the ICM structure factor and abso- 
lute value of P. Within the simulation errors, the two 
order parameters have the same ordering temperature, 
demonstrating the strong mutual influence of ferroelec- 
tricity and magnetic order. The corresponding temper- 
ature phase diagram is presented in Fig. (Sfb). Here, we 
find only one ICM phase with a direct phase transition 
from the disordered state. In order to account for the 
collinear ICM phase at the intermediate temperatures, 
one has to include in the model the anisotropy associated 
with the orthorhombic distortions. Then, different crys- 
tallographic axes become inequivalent and the correla- 
tions of the corresponding spin components become finite 
at different temperatures, due to magnetic anisotropy. 33 

IV. CONCLUSIONS 

We have demonstrated that the DMI provides a nat- 
ural explanation for the coexistence and strong coupling 
between ferroelectricity and IC magnetism in the per- 
ovskite manganites. The empirically estimated value 
of 7 is in agreement with previous experiments. The 
ICM-FE phase appears in the nearest-neighbor interac- 
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FIG. 5: (Color online) Finite temperature MC results for 
Hamiltonian H; ki = 1, k% = 10. (a) Typical temperature 
dependence of the ICM and FE order parameters. The order 
parameters are normalized to their maximum value at T — 0. 
(b) Temperature phase diagram. 



tion model (|10|) as a result of the competition between 
the double exchange, superexchange and DMI. Therefore, 
our model favors a mechanism of stabilization of helical 
magnetic order which is qualitatively different from the 
theories based on spin frustration>i2iiL2£ Unlike in the 
model suggested in Ref. the noncollinear spin state is 
realized in the undoped regime, and thus it is stable with 
respect to phase separation^ 

We believe that, due to its universality, the DMI is 
also relevant in other IC multiferroics. Indeed, the DMI 



is independent of the orbital structure, it only involves 
the interaction between spins induced by the symmetry- 
breaking (in particular FE) ionic displacements. In sup- 
port of this assertion, we note that, despite the diversity 
of crystal structures, P 1 k always for all the IC multifer- 
roics discovered so far, with or without applied magnetic 
fieldfSiiSiSiLiiiii^ which is naturally explained by the 
DMI. 
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